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Abstract 

An efficient method of modeling trapezoidal built-up wing structures is developed by coupling, 
in an indirect way, an Equivalent Plate Analysis (EPA) with Neural Networks (NN). Being 
assumed to behave like a Mindlin-plate, the wing is solved using the Ritz method with Legendre 
polynomials employed as the trial functions. This analysis method can be made more efficient by 
avoiding most of the computational effort spent on calculating contributions to the stiffness and 
mass matrices from each spar and rib. This is accomplished by replacing the wing inner-structure 
with an "equivalent" material that combines to the skin and whose properties are simulated by 
neural networks. The constitutive matrix, which relates the stress vector to the strain vector, and the 
density of the equivalent material are obtained by enforcing mass and stiffness matrix equities with 
regard to the EPA in a least-square sense. Neural networks for the material properties are trained in 
terms of the design variables of the wing structure. Examples show that the present method, which 
can be called an Equivalent Skin Analysis (ESA) of the wing structure, is more efficient than the 
EPA and still fairly good results can be obtained. The present ESA is very promising to be used at 
the early stages of wing structure design. 
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Introduction 


Traditionally, the conceptual design is often carried out using simplified relations that have 
been learned from previous similar designs. This approach can only result in incremental 
advancements in technology because experience-based design makes large step design 
extrapolations too risky. This risk can be reduced significantly if the design is based on physics- 
based high-fidelity models since using these models one can predict the consequences of large 
design extrapolations with a greater confidence. 

Physics-based high-fidelity models, such as those in the Finite Element Analysis (FEA) for 
structures. Computational Fluid Dynamics (CFD) for aerodynamic loads etc., are increasingly 
being used as early as possible in airplane design. This is not surprising because it is estimated that 
about 90% of the cost of a product is committed during the first 10% of the design cycle, and an 
accurate analysis is often crucial to obtaining a good estimate of the manufacturing and life-cycle 
costs of a product. 

For structural analysis, the Finite Element Method (FEM) is widely used because of its 
generality, versatility and reliability, but its use in the early design stages faces some major 
obstacles: A prohibitive preparation time for the FE model data, and a large amount of CPU time 
for problems with a high number of degrees of freedom. This is especially true for complex built- 
up structures such as the airplane wings. 

In view of this situation, Kapania and Liu' have developed an efficient method, using the 
equivalent plate analysis (EPA), for studying the static and vibration problems of general 
trapezoidal built-up wing structures composed of skins, spars, and ribs. The model includes the 
transverse shear effects by treating the built-up wing as a plate following the Reissner-Mindlin 
theory (FSDT). The Ritz method is used with the Legendre polynomials being employed as the trial 
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functions. The Legendre polynomials have the orthogonality property, and this is an advantage over 
using simple-polynomials as the trial functions, which is known to be prone to numerical ill- 
conditioning problems. 

But there is a problem in the method in Kapania and Liu 1 , that is, since complex trial functions 
are used, the stiffness and mass matrices cannot be integrated analytically, as in the case of Livne 2 . 
and numerical integration over every item of the structural components takes a large part of the 
computing effort. For instance, in solving a free vibration problem, evaluating various matrices 
requires about 2/3 of the total CPU time. 

This problem is addressed in this paper. Instead of evaluating the matrices over all components 
of the wing structure, evaluation is performed only over the skins, whose "equivalent" material 
constitutive matrix and mass density distribution are changed accordingly to incorporate the spar 
and rib effects. The new skin material properties are simulated using Neural Networks in terms of 
the wing design variables. As shall be shown, while the Neural-Network-aided EPA, which will be 
called Equivalent Skin Analysis (ESA) hereafter, gives almost equally good results, it uses only a 
fraction of the CPU time spent in the ordinary EPA in evaluating the matrices. 

An Equivalent Plate Analysis (EPA) Method 

The Reissner-Mindlin method, a First-order Shear Deformation Theory (FSDT), is based on 
two assumptions for the displacement field: (1) A straight line normal to the non-deformed middle 
surface remains a straight, but not necessarily normal to the deformed surface, line after 
deformation; and (2) The transverse normal stress can be neglected in the constitutive relations. 

We want to analyze a trapezoidal wing by assuming that it behaves as a plate whose 
deformation satisfies the Reissner-Mindlin displacement field. For the convenience of calculation, a 


3 



coordinate transformation from (x, y) to (%,?!) is performed, which transforms the wing plan 
surface from the original trapezoidal to a square. Details of the coordinate transformation are 
available in Ref. 1 . 

By representing the displacement components as the sum of terms )P , ( rj ) , where P s 
are the Legendre polynomials, i = 1, ■■■,/, j = \,---,J {I,J are integers), we can call the vector 
composed of all the coefficients of the terms as the generalized displacement vector. Based on this 
displacement vector and considering the total strain and kinetic energy of the wing structure, we 
can derive the stiffness and mass matrices as in the following forms: 


[K] = ljj{CY[T) T [D][T}[C]dV 

V 

(1) 

[M] = \l\p[H] T [ZZ)[H]dV 

(2) 


V 


where the integration domain V includes all and only the spaces which the components of the wing 
occupy. Details of obtaining matrix [C] etc. can be found in Ref. 1. For different wing components, 

skin, spar or rib, the material constitutive matrix should be treated differently 1 The integration in 
the (q,r}) plane is performed using the Gaussian quadrature. 

The boundary conditions are approximated by applying springs with very large magnitudes of 
stiffness on the boundaries. For clamped edges both linear and rotational springs are needed, and 
for simply supported edges only linear springs are used. An integration on the boundaries gives 
r K 1 . The total stiffness and mass matrices of the wing structure is given as 

l spring * 

IK ] = ] + {K, n „], \K ] = lK, k „ } + [K„, } + [K„ b ] , ( 3 ) 
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The natural frequencies and mode shapes for the free vibrating wing can be determined by 
applying the Lagrange’s equations, which results in the following eigenvalue problem 

}=0 (5) 

where T = <2T is an eigenvalue of the system of equations, ft) is the corresponding frequency in 
radians/second, and {x} is the corresponding eigenvector. 

Static problems (including stress distributions) can be readily solved by making use of the total 
stiffness matrix [K miul ]. For details one can refer to Ref. 1. 

Compared with the methods in Kapania and Lovejoy 3 and Cortial 4 , the formulation in Ref. 1 is 
such that there is no limitation on the wing thickness distribution. As shown in Ref. 1, the method 
shows a good performance for both static and vibration problems in comparison with the FEA 
using MSC/NASTRAN. 

As Eqs. (1) and (2) are being evaluated by numerical schemes such as the Gauss quadrature 1 , 
the major computing effort is spent in calculating the contributions from the large number of the 
local components, namely, spars and ribs. The purpose of the present paper is to reduce the 
computing effort in obtaining [AT] and [M] by incorporating the contributions from those local 
components to the global components such as the skin. This is accomplished by postulating a new 
material constitutive matrix and mass distribution for the skins and relating the terms of the 
constitutive matrix and mass density to various design variables of the wing structure by using 
Neural Networks, as shall be shown in the following sections. Design variables for wing structure 
can be sizing-type variables (skin thickness, spar or rib sectional area etc.), shape variables (the 
plan surface dimensions and ratios), and topological variables (total spar or rib number, wing 
topology arrangements etc.). 
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Application of Neural Networks in Structural Problems 


Artificial Neural Networks (ANN), or simply Neural Networks (NN) are computational systems 
inspired by the biological brain in their structure, data processing and restoring method, and 
learning ability. More specifically, a neural network is defined as a massively parallel distributed 
processor that has a natural propensity for storing experiential knowledge and making it available 
for future use by resembling the brain in two aspects: (i) Knowledge is acquired by the network 
through a learning process; (ii) Inter-neuron connection strengths known as synaptic weights (or 
simply weights) are used to store the knowledge 5 . 

Major steps of utilizing NN include: (i) specifying the topology or the structural parameters 
(number of layers, number of neurons in each layer, etc.) of the NN, (ii) training of the NN, 
corresponding to the learning process of the brain, (iii) simulation, corresponding to the recalling 
function of the brain. 

The NN has the following properties: (i) Many of its kind are universal approximators, in the 
sense that, given a dimension (number of hidden layers and neurons of each layer) large enough, 
any continuous mapping can be realized; Therefore, (ii) a NN provides a general mechanism for 
building models from data, or gives a general means to set up input-output mapping (iii) The input 
and output relationship of NN can be highly nonlinear; (iv) A NN is parallel in nature, and it can 
make computation fast when executed in a parallel computer, though NN can be simulated 
in ordinary computers in a sequential manner. 

Neural Network has found numerous applications in science and engineering, from biological 
and medical sciences, to information technologies such as artificial intelligence, pattern 
recognition, signal processing and control, and to engineering areas as civil and structural 
engineering. In the field of structural engineering, there have been a lot of attempts and researches 
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making use of NN to improve efficiency or to capture relations in complex analysis or design 
problems. The following are a few examples. Hajela and Berke 6 gave an overview of the neural 
computing approach, and examined the role of neural computing strategies in structural analysis 

and design. Abdalla and Stavroulakis 7 applied NN to represent experimental data to model the 
behavior of semi-rigid steel structure connections, which are related to some highly nonlinear 
effects such as local plastification etc. Several cases of neural network application in structural 

g 

engineering can be found in Vanluchene and Sun . All the problems treated in Ref. 8 had been 

reproduced in Gunaratnam and Gero 9 with a conclusion that representational change of a problem 
based on dimensional analysis and domain knowledge can improve the performance of the 

networks. In Liu, Kapania and VanLandingham 10 , methodologies of applying Neural Networks 
and Genetic Algorithms to simulate and synthesize substructures were explored in the solution of 1- 
D and 2-D beam problems. 

For the efficient simulation of the structural performances of complex wings, there can be two 
directions to apply NN as specified in the following: 

(1) Direct Application 

In this case, the input layer includes all the design variables of interest (for instance, the four 
shape parameters of the wing plan form). The output layer gives the desired structural 
responses, such as natural frequencies etc. The EPA is being used as the training data generator, 
though if necessary, results obtained using the FEM can also be used as the training data. 
Preparation of training data is very important, and the training algorithm used also greatly 

impacts the process of training 1 1 . Caution must be taken in specifying the network parameters 
and training criterion, such that the results of the trained network would not oscillate around the 
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training data. Once the networks are trained, structural responses at any design point can be 
recalled in a fraction of a second and this is really favorable in a design situation 12 . 

(2) Indirect Application 

Here it is desired to find a way of incorporating NN into the application of the equivalent plate 
analysis (EPA) of complex wing structures, other than just making use of results generated by 
EPA as the training data base. Note that in the EPA of a complex wing, the computation effort 
is mainly spent on integrals for generating the contribution from the inner-structural 
components of the wing, the spars and the ribs, to the stiffness and mass matrices. If an 
equivalent anisotropic material can be found to replace the inner components, in terms of an 
equivalent skin, such that the new composite wing has similar global properties as the original 
one, then the EPA can be performed more efficiently. Solution of the adequate material 
properties of the anisotropic material is the major obstacle here. The role of NN will be relating 
the material properties to all kinds of wing design parameters, when there exists enough data 
base for training. This way of applying NN has been claimed to be the best use of the Neural 

9 

Networks in structural engineering : 

"The real usefulness of neural networks in structural engineering is not in reproducing existing 
algorithmic approaches for predicting structural responses, as a computationally efficient 
alternative, but in providing concise relationships that capture previous design and analysis 
experiences that are useful for making design decisions. " 

This is the path that we followed in this study. 
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Equivalent Skin Analysis (ESA) of Wing Structures 

Since the calculation of various integrals in Eqs. (1) and (2) are time-consuming, it is desired to 
replace the actual wing structure by an equivalent continuum model, that is, one that is composed 
of a skin-like material, whose constitutive matrix [D] and distribution of mass in Eqs. (1) and (2) 
are to be decided. 

It is assumed that the mass density p is a function of position in the plan form while each term 
of [D] is a constant throughout the wing area. There can be other choices, as will be discussed 
later. 

We are going to solve the above problem by requiring that the stiffness and mass matrices of 
the equivalent model are most approximate to those of the actual wing in a least squares sense. This 
gives the following procedures. 

(a) The Constitutive matrix 

Let’s write [ K ] = [K slrain } = [K t) ] as the target matrix, and the stiffness matrix of the equivalent 
continuum model is 

[k] = ^g m gAC] T mn [711 [0][n m „ [C] mn 

" _ ( 6 ) 

= [k :l ] = [^D M (^Gi mn )] = ^DJG; q ] 

p q m n p,q 

where g m and g n are the Gauss quadrature weights; the constitutive matrix [D] relates the stress 
and strain vectors by {a } = [D]{f } , and D pq is the p -th row, 4-th column term of the constitutive 
matrix; 

m ( m = \, w ",K) corresponds to the m -th Gauss integration position in the x-direction. 

K is an integer with a usual value of 6 or 8; 

n ( n = \,-,K) corresponds to the n -th Gauss integration position in the y-direction; 
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p — 1, ■ ■ ' ,5 is the row number of [D] ; 
q = 1, - ,5 is the column number of [D]; 


G^ mn =g m g„[Cy mn [T] 


T rT'i?' 

mn 


0 ••• l 

0 

1 q 


0 


P [T] mn [C] mn ,G; q =^G^ n 

• m,n 

5 


By constructing an error function 

E, =1 


'•y 


where w* are the weight coefficients, and by requiring 


dE 




3D 


pq <■! 


dD 


p<i 


and noting 


dK„ 


dD 


= G' 1 . ( G'j, q means the i, j -th term of matrix [G pq ] ), 


P‘7 


(7) 


( 8 ) 


we can obtain the constitutive matrix term [D pq ] by solving the following linear equation set: 


Y(Y W K G" G '' )D = Y W H K„G 
q pq pq 7 pq v v 


•j 

pq 


(9) 


p q ij 




i,j = \, -,K and p,q,p\q'= 1, ■ ■ • ,5 . 

This is an equation set with 25 unknowns. Since i, ; = 1, • • • , N and N = 5K 2 is usually very 
large (if we use the Legendre polynomials of 6 terms as the basis functions, N = 180 , if 8 terms are 
used, then N = 320), the job of generating the matrix in Eq. (9) will be very extensive. 
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If [D] is assumed to be symmetrical, then Eq. (9) will become 


IS 


VV, 




,.)G" ]G" 

p q q p pq 


D ,,=I 


' <J 




( 10 ) 


<r ‘i M 


where (p\q’) and ( p,q ) have the following 15 instead of 25 combinations: (1,1), (1,2), (1,3), 
(1,4), (1,5), (2,2), (2,3), (2,4), (2,5), (3,3), (3,4), (3,5), (4,4), (4,5), and (5,5). 


(b) Mass distribution 

Let’s write [ M ] = [M mml ] = [M tj ] as the target matrix, and 


[M]=IIP„.«,^[ff]L[z]L[zUff]. 1 .=IIP™[f]., on 

m n m n 

as the mass matrix of the continuum model, where m (m - 1, • - • , K) corresponds to the m -th Gauss 
integration position in the x-direction, rc ( /t = 1, • • • , K - ) corresponds to the n-th Gauss integration 
position in the y-direction, p mn is the mass density of the equivalent model at position ( m,n ), g m 
and g n are integration weights, and =[H] T mn [Z] T mn [Z] mn [H] mn is a NxN matrix varying 

with position ( m.n ). 

By constructing an error function 

E„ =S<[w, ,(/>_)-«, T (12) 

where vv y w are weight coefficients, and by requiring 


dpm. 




•-J 


dM: 

Jj 

Zp m 


= 0 , 


and noting 



( is the i, j -th term of matrix F mn ), 
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we can obtain the mass distribution p mn by solving the following linear equation set with K : 
unknowns: 

X ( X w >! F mn F L )P m w = X w " m',n'= 1, — » K . (13) 

m .n ' i,j i.j 

In the present study the following weight coefficients are used: 

w,f = 1 0| ATf " I /max ( I " I) + 1 

i , / ' 7 i , ( 14 ) 

wf = 1 o| M f " | / max ( | A/ f m | ) + 1 

The basic idea behind this choice is that we want to form the equivalent matrices more in the way 
of the skin’s, which is more like a plate than the other components of the wing, i.e. spars and ribs. 

Several choices about the variation of p and [D] have been tried, but it is found that the 
present assumptions give the best results in terms of feasibility and accuracy. For instance, to be 
consistent with the assumption that each term of [D] is a constant throughout the wing area, p can 
also be assumed a constant. This certainly decreases the accuracy of the method due to the loss of 
flexibility in varying p to simulate the target mass matrix [M], but the resultant reduction in 
computational effort is small since in the first place, forming Eq. (13) and training the p -related 
neural networks do not need much CPU time. In other cases, [D] was assumed to be variable in the 
span-wise direction or throughout the wing area, but it is found that although the equivalent 
material is more flexible to simulate the target stiffness matrix [AT] , the resultant [K] usually has a 
larger abstract error and the solution of the free vibration problem usually gives worse natural 
frequencies. Moreover, the CPU time needed for generating Eq. (10), which requires the major 
computational effort in our method, increases in a factor of about K (number of Gauss integration 
points, usually with a value of 6) in the case of [D] being variable in the span-wise direction. In the 
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case of [ D ] being variable throughout the wing area, the increase can be as large as K : times. As 
we shall see in the following examples, these kinds of increase in CPU time are formidable. 

Examples and Discussion 

The Neural-Network-aided equivalent plate analysis, or briefly Equivalent Skin Analysis (ESA) 
method is compared with the ordinary equivalent plate analysis (EPA) 1 in three cases where 4 to 6 
design variables are involved respectively. In some of the results, FEA results employing 
MSC/NASTRAN are also provided as a benchmark. 

Some common parameters of the built-up wing structures will be specified herein if not 
specified otherwise in the following cases. The sections were generated by the Karman-Trefftz 

transformation 13 . The section thickness-chord ratio varies from 0.15 at the root to 0.06 at the tip. 
Skin Thickness t 0 =. 11 8 in ; Spar Cap Height h x = .197m , Spar Cap Width /, = .373 in , Spar Web 

Thickness r, = ,058m (for definition of h x etc. one can see Ref. 1 ); the ribs have the same cap 
dimensions and web thickness as the spars. Mass Density p - 2.526 x 10 lb ■ sec"/ in'* , Young’s 

Modulus E = 1.025x 10 1 lb/in 2 , Poisson's Ratio v =0.3. The wing is clamped at the root. 

(a) Four-variable case: design space I 

In this case 4 spars and 10 ribs are evenly distributed inside the wing plan form under the skins. 
For a trapezoidal wing, there are four major independent shape variables: sweep angle A , aspect 
ratio a , taper ratio t , and plan area A (see Fig. 1). A 3 4 full factorial experimental design with 3 
levels in A , Ct , I and A respectively, was used. Particulars of the levels of every variable are as 
follows: 

A = {O’ ,15’ ,30" } , a = {1.0,1.75,2.5}, j 
r = {0.3,0.45,0.6}, A = {2000,3500,5000}m'.J 
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For each point in this design space, the EPA is carried out, then Eqs. (10) and (13) are used to 
generate the 15 constitutive matrix terms and mass densities at 36 (6x6) Gauss sampling points, 
and the ESA is performed. For each of these parameters, a feed-forward neural network with a 
structure of 4x15x10x1, i.e. 4 inputs, 15 neurons in the first hidden layer, 10 neurons in the 
second hidden layer, and 1 output, is trained using the MATLAB NN Toolbox function trainlm that 
trains feed-forward network with the Levenberg-Marquardt algorithm 10 . Therefore, there are totally 
15+36=5 1 networks to be trained. There are 8 1 ( 3 4 ) sets of training data, which are non- 
dimensionalized before the training process. Once the networks are trained, the input-output 
relationships can be readily retrieved by using the MATLAB function simuff. 

The major computational effort was spent in generating the 81 sets of training data, with about 
45 hours of CPU time being spent on a PII/350 personal computer, while less than 1 hour of CPU 
time being used in training the neural networks. A set of results are given in Fig. 2 where 25 points, 
which mean 25 new designs, were randomly chosen within the design space box. Upon each new 
design both the EPA and the ESA are performed. The plan forms of the new design are shown in 
Fig. 2(a). The first 10 natural frequencies by the EPA and the ESA are compared in Fig. 2(b) and 
their relative differences (based on the EPA results) are shown in Fig. 2(c). It can be seen that 
except for a very few cases (3 out of 250), the relative difference is within - 10%~ 10%. 

Fig. 3(a) shows 16 new designs through a randomly chosen path inside the design space box 
which is defined as 


v' = v'(l - a 1 ) + v(a' ,j = 1, - • • ,4 
v 1 = A, v 2 = or, v 3 = r,i / 4 = A 
a' =s n ' ,n i - r j /(I - r ] ). 


(15) 
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where v 0 ' and v( are the lower and upper bounds of variable v' , for instance, v ] 0 =0°. vj = 30’ 
etc., se [0.1] is the range of a shape variable, and r J (j = l,--4) are randomly determined values 

between 0 and 1. Results of natural frequencies of the first 6 modes for wing structures defined by 
points along a path with n , =0.945, n 2 =8.200, n 3 =3.203 and n 4 =1.778 are shown in Fig. 3(b), 
where it can be seen that results by the EPA and the ESA agree with each other quite well. 

While Figs. 2 and 3 are about free vibration frequencies, Fig. (4) shows some static results. For 
an arbitrary new design whose plan form as shown in Fig. 7(a), a down-ward (- 2 direction) point 
force of 1 lb is applied at the mid-point of the wing tip (actually the force is divided into 
components acting on the two spar tips close to the mid-point). Fig. 4(a) shows displacements 
along the leading-edge by the EPA and the ESA, where u,v,w are displacement components in the 
chord-wise, span-wise, and vertical directions respectively. Figs. 4(b) and 4(c) show the Von Mises 
stress distributions at the wing root and the central span-wise line respectively. It can also be seen 
that the EPA and the ESA give very comparable results for the static case. 

(b) Six-variable case: design space II 

In this case spars and ribs are evenly distributed inside the wing plan form but their numbers are 
design variables. A 3 6 full factorial experimental design with 3 levels in A , a , t , A and 
numbers of spars and ribs, n !par and n rib respectively, was used. Particulars of the levels of every 

variable are as follows: 

A = {0,15,30}, or = {1.0,1.75,2.5}, 
r = {0.3,0.45,0.6} , A = {2000,3500,5000}w 2 , > 
n var = {2,4,6}, n nh = {7,10,13}. 

Similar to case I, for each point in the design space II, the EPA is carried out and Eqs. (10) and 
(13) are used to generate the 15 constitutive matrix terms and the 36 mass densities which are then 
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used to perform the ESA. A total of 5 ] feed-forward neural networks with a structure of 
6x15x10x1 are trained using the MATLAB NN Toolbox function trainlm. There are 729 ( ) sets 

of data that could be used for training, but it was found that at some design points the differences 
between the natural frequencies by the EPA and the ESA become too large. Therefore, a screening 
process was introduced, in which any point where the maximum relative difference between the 
first 10 natural frequencies by the EPA and the ESA surpasses 20% will be discarded. This led to 
the removal of 28 points through the process, therefore 701 sets of data were used for training. 

Generating the 729 sets of pre-training data used about 152 hours of CPU time on the Crunch 
(SGI Origin 2000 with eight R 10000 processors) of the College of Engineering, Virginia Tech, and 
training the neural networks spent about 2 hours on a PII/350 PC. A set of results are given in Fig. 

5 where 25 points were randomly chosen within the design space box. The plan forms of the new 
designs are shown in Fig. 5(a), where dashed lines indicate the spar or rib positions. The first 10 
natural frequencies by the EPA and the ESA are compared in Fig. 5(b) and their relative differences 
(based on the EPA results) are shown in Fig. 5(c). It can be seen that the relative difference is 
within -5%~15%. 

Fig. 6(a) shows 16 new designs through a randomly chosen path inside the design space box 
which is defined as 


v 1 =v£(\-a') + v(a i ,j = \,---,6 

v 1 = A, v 2 =a,v } = r,v 4 = A, v 5 =n vpar ,v 6 =n nh ■ 

a 1 = s n ' ,n } = r J /(l- r y ). 


(16) 


where r (j = l, ••■,6) are randomly determined values between 0 and 1, and see Eq. (15) for the 
definition of other symbols. Results of natural frequencies of the First 6 modes for wing structures 
defined by points along a path with n , =0.2243, n , =0.8591, n , =0.2064, n 4 =3.0700, n , =2.21 96 


16 


and n 6 =0.9440 are shown in Fig. 6(b), where it can be seen that results by the EPA and the ESA 
agree with each other quite well. 

Now some static results. For an arbitrary new design whose plan form is shown in Fig. 7(a), a 
down-ward (- z direction) point force of 1 lb is applied at the mid-point of the wing tip. Fig. 7(b) 
shows displacement components along the leading-edge by the EPA and the ESA, compared FEA 
using MSC/NASTRAN. Figs. 7(c) and 7(d) show the Von Mises stress distributions at the w'ing 
root and the central span-wise line respectively together with the FEA results. Comparison of the 
natural frequencies of this wing as given by the EPA, the ESA and the FEA is shown in Table 1 . It 
can be seen that the EPA and the ESA results are close, and they all agree quite well with the FEA 
results. 

(c) Design space III 

In this case a wing plan with A = 30° , s = 192 in,b = 12 in, and a = 36 in (see Fig. 1 for 

definitions of 5 , b , and a ) is used. A 2 4 x 3 2 full factorial experimental design with 2 levels in t 0l 

(skin thickness at wing tip), a n (skin thickness increment ratio at root over the tip), h x (spar cap 

height) and h 2 (rib cap height), and 3 levels in n spar and n rih , is carried out. The skins are assumed 

to vary linearly from the root to the tip. Particulars of design space III are as follows: 

t 0l = (1,3} x 0.1 18 in, a n ={ 0,2}, 

h t ={1,3} x 0.197 in, h 2 = {1,3} x 0.197 in, > 
n spar = (2, 4,6), n nh = (6,10,14}. 

There are 144 sets of data for training. Generating these data sets used much less CPU time than 
in the case of design space II (about 30 hours). A set of results are given in Fig. 8 where 16 points 
were randomly chosen within the design space box. The plan forms of the new designs are shown 
in Fig. 8(a), where dashed lines indicate the spar or rib positions, and the skin thickness at the wing 
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root and tip. and cap heights of spars and ribs are represented as shown in Fig. 10(a). The first 10 
natural frequencies by the EPA and the ESA are compared in Fig. 8(b) and their relative differences 
(based on the EPA results) are shown in Fig. 8(c). It can be seen that the relative difference is 
within -5 c 7c~\5 c /c. 

Fig. 9(a) shows 16 new designs through a randomly chosen path inside the design space box 
which is defined as 


v ' = Vq ( 1 - a 1 ) + v(a ' , j = 1, • • • ,6 

v' =t 0l ,v 2 =a n ,v- =/t,,v 4 =h 2 ,v 5 =n spar ,v 6 =n rih > 

a' = s"’ ,n j = r /(l - r y ). 


(17) 


where r (j = 1, • ,6) are randomly determined values between 0 and 1, and see Eq. (15) for the 
definition of other symbols. Results of natural frequencies of the first 6 modes for wing structures 
defined by points along the path with n , =0.003 1 , n 2 =0.9999, n 3 =0.2089, n 4 =64.7024, n 5 =0.9067 


and n 6 =0.5325 are shown in Fig. 9(b), where it can be seen that results by the EPA and the ESA 
agree with each other quite well. 

For an arbitrary new design whose plan form is shown in Fig. 10(a), a down-ward (- z 
direction) point force of 1 lb is applied at the mid-point of the wing tip. Fig. 10(b) shows 
displacement components along the leading-edge by the EPA and the ESA, compared with FEA 
using MSC/NASTRAN. Figs. 10(c) and 10(d) show the Von Mises stress distributions at the wing 
root and the central span-wise line respectively together with the FEA results. Comparison of the 
natural frequencies of this wing as given by the EPA, the ESA and the FEA is shown in Table 2. 
A^ain, it can be seen that the EPA and the ESA results are close, and they all agree quite well with 
the FEA results. It is noted that a coarser design space III does not worsen the accuracy of the ESA. 
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The CPU time savings using the ESA are obvious. For instance, when 6 terms of the Legendre 
polynomials {K - 6 ) are used, about 85% less CPU time is spent in evaluating the stiffness and 
mass matrices compared with the EPA, where matrix evaluation takes about 68% of the total CPU 
time when solving the free vibration problem. When K = 8, about 83% less CPU time is spent in 
evaluating the matrices compared with the EPA, where matrix evaluation takes about 65% of the 
total CPU time. Generally speaking, the results given by the ESA in design space II and HI are as 
good as those in design space I although the number of variables increases from 4 to 6. 

Conclusion 

An efficient method for the analysis of built-up wing structures using Equivalent Skin Analysis 
(ESA) has been developed. The CPU time spent on evaluating the matrices, which takes about 
63-68% the total CPU time in the EPA applying to solve a free vibration problem, is about 5-6 
times less when the ESA is used. Three groups of examples with 4 or 6 design variables show that 
fairly good results can be obtained. This method is very promising to be used at the early stages of 
wing design. 
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Table 1 Natural frequencies (rad/sec) of the wing in Fig. 7(a) 


Mode No. 

1 

2 

3 

4 

5 

Mode Shape 

1 st bending 

2 nd bending 

1 st torsion 

In plane 

2 nd torsion 

EPA 

279.3 

982.8 

1057.5 

1447.4 

1945.5 

ESA 

274.5 

984.1 

1045.1 

1440.3 

1936.3 

FEM 

279.9 

965.6 

973.5 

1454.4 

1830.8 


Table 2 Natural frequencies (rad/sec) of the wing in Fig. 10(a) 


Mode No. 

1 

2 

3 

HHHHMi 

5 

Mode Shape 

I s ' bending 

■UBB3I 



3 rd bending 

EPA 

71.9 

233.9 


452.2 

479.9 1 

ESA 

70.8 

239.4 

358.4 

469.4 

504.8 : 

FEM 

66.0 

222.6 


413.1 

468.0 1 








Fig. 1 Plan configuration of a trapezoidal wing: A = -js(a + b),a = s 2 j A,r = ajb. 

PDP 

PPD^ 



Fig. 2 (a) 25 randomly chosen wing plan forms in design space I 
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Fig. 3 (a) 16 wing plan forms systematically varying through design space I 



Fig. 3 (b) Comparison of the first 6 frequencies by EPA and ESA 




O.OE+OOI 


£-1 OE-04 

O 

c 


> 

5 -2. 0E-04 



0 1 0 20 30 40 50 60 70 80 

Distance from root (inch) 

Fig. 4 (a) Comparison of displacements by EPA and ESA under 1 lb tip force 
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Fig. 4 (b) Comparison of the Von Mises stress at wing root by EPA and ESA under 1 lb tip 





Von Mises stress (psi) 



Fig. 4 (c) Comparison of the Von Mises stress along central span-wise line by EPA and ESA 

under 1 lb tip force 



Fig. 5 (a) 25 randomly chosen wing plan forms in design space II 
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Fig. 6 (a) 16 wing plan forms systematically varying through design space II 



Fig. 6 (b) Comparison of the first 6 frequencies by EPA and ESA 
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Fig. 7 (c) Comparison of the Von Mises stress at wing root by EPA and ESA under 1 lb tip 

force 



Fig. 7 (d) Comparison of the Von Mises stress along central span-wise line by EPA and ES A 

under 1 lb tip force 
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Fig. 9 (b) Comparison of the first 6 frequencies by EPA and ESA 



Fig. 10 (a) An arbitrarily chosen wing design in design space III 
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Fig. 10 (b) Comparison of displacements by EPA and ESA under 1 lb tip force 



Fig. 10 (c) Comparison of the Von Mises stress at wing root by EPA and ESA under 1 lb tip 

force 
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Fig. 10 (d) Comparison of the Von Mises stress along central span-wise line by EPA and ESA 

under 1 lb tip force 
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